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SECTION  1 


INTRODUCTION 

One  of  the  primary  considerations  in  the  design  of  structural  sandwich 
components  is  the  stability  of  the  panel  when  subjected  to  inplane  loadings. 
These  applied  loads  consist  of  edgewise  compression,  edgewise  shear  and 
edgewise  bending  moments,  as  illustrated  in  Figure  1. 

The  amount  of  literature  available  concerning  the  buckling  of  flat, 

retangular  sandwich  panels  is  rather  limited.  There  are  several  references 

1 # 

to  uniaxial  edgewise  compression  and  several  references  to  the  problem 

2 

of  shear  alone  . The  case  of  combined  edgewise  bending  and  compression 
has  not  been  extensively  investigated.  Other  observations  have  considered 
the  effects  of  compression  combined  with  shear  by  the  use  of  interaction 
formulas  . The  results  presented,  therefore,  have  very  little  in  the  way 
of  previous  results  that  would  in  some  way  support  the  material  presented. 

This  report  presents  a method  for  the  prediction  of  the  general 
stability  of  flat,  clamped  sandwich  panels  loaded  by  combined,  edgewise 
bending,  biaxial  compression,  edgewise  shear,  and  uniaxial  compression. 

It  is  anticipated  that  the  method  presented  will  be  a useful  tool  in  the  design 
and  analysis  of  sandwich  components. 

1.  1 DEFINITION  AND  SCOPE 

The  analysis  presented  considers  buckling  of  flat,  rectangular 
sandwich  panels  subjected  to  combinations  of  uniaxial  edgewise  compression, 
biaxial  edgewise  compression,  edgewise  bending,  and  edgewise  shear  loading. 
The  component  layers  are  restricted  to  linear  elastic  behavior  and  small 
displacements  for  the  clamped  panels  being  considered. 


>!< 


Numerical  superscripts  indicate  References. 


COMBINED  LOADS 
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The  core  is  considered  to  be  completely  rigid  in  the  direction  normal 
to  the  panel.  The  resistance  of  the  core  to  extension  and  bending  is  assumed 
to  be  negligible;  therefore,  the  strain  energy  may  be  taken  to  consist  only  of 
contributions  due  to  transverse  shear  moduli  of  the  core  which  may  be  helpful. 

The  face  sheets  are  assumed  to  be  isotropic  but  may  have  different 
properties  and  thicknesses.  Critical  combined  loads  of  classical  plates 
may  be  obtained  as  a special  case  (Love -Kir choff  hypothesis),  since  the 
flexural  rigidity  of  the  faces  is  not  neglected. 

1.  2 METHOD  OF  SOLUTION 

Linear  analyses  of  sandwich  deformations,  in  general,  fall  into 
4 

two  categories  as  follows:  formulations  in  terms  of  face  displacement, 

and  analyses  in  terms  of  functions  related  to  core  displacements. 

Theoretically,  the  two  methods  are  equivalent,  but  with  regard  to  the 
boundary  conditions,  permit  slightly  different  numerical  treatments. 

3 

The  results  presented  are  based  on  the  so-called  "tilting  method"  , 
by  which  the  displacements  in  each  face  are  related  to  the  displacement  of 

!the  core.  The  solution  is  obtained  by  using  the  Ritz  method  to  minimize 

the  total  potential  energy  in  the  panel.  For  the  purpose  of  uncoupling  the 
B 

flexural  and  extensional  strain  energies,  it  is  assumed  that  there  exists  a 

. 

"neutral  plane"  in  which  no  inplane  displacements  occur  during  the  buckling 

i 

mode  deflection.  The  neutral  plane  location  which  corresponds  to  the 
given  buckling  mode  is  then  treated  as  an  additional  degree  of  freedom 
in  the  minimization  of  the  potential  energy. 

1 

It  is  clearly  noted  that  a linear  analysis  is  only  valid  under  certain 
ideal  conditions.  The  most  important  conditions  are  that  prior  to  buckling, 
the  structural  components  experience  only  infinitesimal  changes  in  geometry, 
and  that  there  is  sufficient  symmetry  such  that  the  assumptions  aimed  at 
uncoupling  modes  of  deformation  are  justified  (i.  e.  , the  existence  of  neutral 
planes).  Also  sandwich  panels  which  have  severely  unbalanced  face  construction 
are  likely  to  require  a more  rigorous,  nonlinear  treatment,  since  the 
assumptions  concerning  the  neutral  planes  may  be  exceeded. 
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SECTION  2 


THEORETICAL  APPROACH 


To  construct  the  total  potential  energy  of  the  panel,  a displacement 
formulation  is  used.  The  theory  on  which  this  analysis  is  based  is  presented 
in  Reference  4.  For  completeness  the  essential  equations  are  presented  here. 

2.  1 GEOMETRY 

In  Figure  2,  the  panel  configuration,  reference  axes  and  loading  con- 
ventions are  indicated.  The  x-y  plane  is  chosen  to  coincide  with  the  panel's 
midplane.  The  stress  resultants  and  displacements  are  shown  to  be  of 
positive  sense  in  the  figure. 

In  Figure  3,  the  assumed  geometry  of  deformation  is  shown  in  the 
xz  coordinate  planes.  A transverse  deflection  occurs,  accompanied  by  trans- 
verse shear  deformations  in  the  core  as  the  panel  buckles.  During  this 
deflection  it  is  assumed  that  there  is  some  plane  located  at  z - e^  where  there 
is  no  further  inplane  displacement  and  that  the  rotations  of  the  panel  in  the  xz 
plane  occur  about  points  in  this  plane  whereas  rotations  occur  around  the  plane 
z = e in  the  yz  plane.  With  these  assumptions,  each  of  the  inplane  displace - 

y 

ments  of  the  plane  can  be  expressed  in  terms  of  the  transverse  deflection  w 
and  the  core  rotations  0 and  Ip  in  the  xz  and  yz  planes  respectively. 


2.2  POTENTIAL  ENERGY 

The  assumed  geometry  (Figure  3)  permits  the  rotations  and  displacement 
about  the  neutral  plane  to  have  the  following  forms  (see  Reference  4 for  details): 


and 


0 = 0 (x,  y).  Ip  = >P  (x,y),  w = w (x,  y) 


3 = . - ft 
Pf  2 f 

V = L +<-')' 


a = P-  +(-l)  e ; f = 1,2 
fy  f y 


(1) 


4 


/ 

FACE  2 (ISOTROPIC) 
E-Z  iv2  $2  ~ != 2 

2(l  + i'2) 


FACE  l (ISOTROPIC) 
E|  G|  = E| 

4 2 ( I +pi) 

CORE  (ORTHOTROPIC) 

Gcxz  » ^cyz 


Figure  2.  Panel  Geometry. 


ure  3.  Geometry  of  Deformation  in  x-z  Plane. 


It  is  again  important  to  note  that  the  neutral  plane  locations  e , e are 

x y 

unknown  parameters  which  will  be  determined  in  the  process  of  minimizing 
the  potential  energy. 

4 

Thus  making  these  substitutions  , the  potential  energy  functional  is 
obtained  as  follows: 


1 | ff  r b r a . 2 2 Uvf  2 

2 f=l  X.  Jo  o fx  f(^’x  2 0,y  } 


0-,  t [ 2 V 0,  + (1  - v )0,  0, 

iy  f 1 x y f v x 


fx  fy 


y x 


2 2 * ~Vf  2 

+ a f,r  t ( 0,  + — 0,  ) 

fy  f y 2 x 


a,  [0.  w»  + VX6,  W,  + ( 1 - V ) 0,  w, 

fx  f r x xx  r x yy  f y 


y xy 


+ a.  t [ 0,  w,  + v,  w,  + ( 1 - vf)  0,  w,  ] 
fy  f y yy  f y xx  f x xy 


1 3 r 2 2 . 2,  , 

+ T t [w,  + 2v,  w,  w,  +w,  +2  (1  - v.)  w,  J ! dxdy 

j t xx  i xx  yy  yy  t xv 


1 r b p a 


2 o o cxz  c 


[ G t (0  + 20  w + w,  ) + G t (b‘  +2^' w,  +w,  )]dxd 

rV7  r 1 v v ' ew-T  « . T-  ir  J 


x x cyz  c 


y y 


+ T J b J a [N  W,  “■  + 2N  W,  w,  +N  w,  z]  dxdy 
2 "o  Jo  x x xy  x y y y J 7 


It  is  important  to  realize  that  the  resultants  N , N , N need  not  be 

x xy  y 

constants  throughout  the  panel,  although  they  must  satisfy  the  equilibrium 
conditions  before  buckling  as  follows: 


7 


N + N =0 

x ,x  xy,  y 

(3) 

"N  + N =0 

*y,x  y.y 


Hence  the  quantities  N and  N may  be  taken  to  be  both  uniform  and 

x y 

linearly  varying  compressive  forces;  having  the  following  forms: 


N 

x 


N 

xo 


+ 


N 

y 


N 

yo 


+ N (x) 
yB 


(4) 


thus  permitting  the  consideration  of  edgewise  bending  forces.  The  special 

case  of  a pure  edgewise  bending  moment  (Figure  1)  will  be  indicated  by  the 

symbols  N and  N such  that 
xB  yB 


NxB<y) 

= NxB 

(1 

- 2 

N (x) 

= N 

(1 

- 2 

yB  ' ’ 

yB 

JL 

b 

X 

a 


(5) 


The  linearly  varying  component  of  the  pure  edgewise  bending  load,  therefore, 
is  twice  the  magnitude  of  the  uniform  component,  and  opposite  in  sign. 
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SECTION  3 


NUMERICAL  SOLUTION 

The  principle  of  minimum  potential  energy  can  be  applied  to  produce 
the  desired  conditions  for  equilibrium.  The  potential  energy  should  be 
expressed  in  terms  of  a complete,  admissible  set  of  functions  that  satisfy 
the  geometry.  Application  of  the  principle  of  minimum  potential  energy 
prodvices  the  appropriate  minimizing  relationships  in  the  form  of  an  eigenvalue 
problem.  By  considering  successively  larger  solution  sets,  each  of  which 
contains  the  previous  one,  monotonic  convergence  to  the  true  solution  is 
surely  guaranteed. 


3.  1 POTENTIAL  ENERGY  BY  THE  RITZ  METHOD 

The  potential  energy  associated  with  the  buckling  of  the  panel  is  given 
by  Equation  2.  This  Equation  is  in  terms  of  w,  0,  and  , thus  admissible 
assumed  mode  functions  must  satisfy  the  conditions  of  continuity  and 

4 

differentiability,  as  well  as  tire  following  imposed  boundary  conditions  for 
tlie  clamped  panel: 


and 


Mo,  y)  = w (a,  y)  = w (x,  o)  = w(x,  b)  = 0 


w,  (o,  y ) = w,  (a,  y)  = w,  (x,  o)  = w,  (x,b)  = 0. 

x x y y 


Therefore,  an  appropriate  set  of  functions  is  as  follow  s: 
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Substitution  of  these  assumed  forms  for  the  displacements  into  the  potential 
energy  (Equation  2)  yields  a function  tt  which  is  thus  expressed  entirely  in 
terms  of  trigonometric  functions  as  in  Equation  6 (see  Appendix  A).  The 
integrations  required  over  the  area  involved  can  be  carried  out  directly, 
according  to  the  following: 


. m-rrx  . rirx  , a 6 

sin sin  dx  = — mr 

a a 2 

m = 0 

rrvrry  riry  , bo  or  b for  _ 

cos — ; — 1 cos  ,L  dy  = — mr  r - U 

b b 2 


. mirx  r-rrx  , 

sin  cos dx  = 

a a 


mux  rirx 

x cos cos 


2 am 


/ 2 2i 

ir(m  - r ) 


mr 

2 


(V] 


„ a r 1 1 

— 6 - ~ [ J + 

4 mr  2 .2 

tt  (m-r)  (m+r) 


2 ^ Amr 


a m=  0 

or  — for 

2 r = 0 


where  6 is  the  Kronec.ker  delta,  and  A is  defined  as  follows: 


mr 


mr 


f 0:  m + r even 
A = ' 

mr  1;  m + r odd 


The  expression  obtained  upon  substituting  the  assumed  modes  for  the 
potential  energy  (Equation  6 ) and  thus  performing  the  indicated  integrations 
as  given  above  is  as  shown  in  Appendix  B. 
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' 


It  is  again  important  to  recognize  that  the  parameters  and  r^. 
contain  the  locations  of  the  neutral  planes  (Equation  1).  Thus  they 
also  are  undetermined  parameters  dependent  upon  the  particular  buckling 
mode;  that  is  , 


0.  , = (a  ) 

fx  fx  mn 


a,  = (cl,  ) 

fy  fy  mn 


Since  the  energy  expression  contains  cubic  terms  such  as  a,  0 which 

fxrnn 

are  evidently  not  quadratic  forms  in  the  unknowns  A , B , C , e , 

mn  mn  mn  xmn 

and  e , the  minimization  of  the  energy  will  involve  a system  of  nonlinear 
ymn 

simultaneous  equations  in  its  present  form.  In  order  to  avoid  this  difficulty, 
it  is  advantageous  to  define  the  following  additional  parameters: 


H = e B 

mn  xmn  mn 


K = e C 
mn  ymn  mn 


Upon  making  this  addition, the  following  parameters  take  on  the  following 


forms  (Equation  1 ). 


, H 

a , . , . f mn 
a fxrnn  @f  + } B 

mn 


, K 

q , , , , i mn 

a.  = P,  + (-D  7= 

fymn  f C 

mn 


? 2 f run 

a.  = a,  a . = + P,(-l)  B 

fxrnn  fxrnn  fx*  p f t mn 


H II. 

mn  *p 

f B B. 

mn  *p 


= ci  f ct  i - ®,  + 3 ( - 1 ) 

fymn  fymn  fy*p  f t 


- K , K. 

' + e,(-n 


K K, 


^mn  ^4p 

, r H 

a c.  B f run 

afxmn  0 fy*p  f f B _ 


+ 3,  (-1) 


f Kl 
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Upon  making  these  substitutions  into  the  expression  for  the  potential  energy 

in  Appendix  B,  the  expanded  and  simplified  expression  for  the  potential 

energy  is  found  in  Appendix  C in  which  all  terms  are  quadratic  in  the  unknowns 

A , B , C , H , and  K 
mn  mn  mn  mn  mn 

A linear  system  of  equations  of  the  unknown  parameters  can  be  formed 
when  the  principle  of  minimum  potential  energy  is  applied. 

3.2  PRINCIPLE  OF  MINIMUM  POTENTIAL  ENERGY 

The  principle  of  minimum  potential  energy  in  the  variational  form  is 
as  follows: 


= o m 

p 

where  tt  is  given  in  Appendix  C.  After  discretization  using  the  Ritz  method, 
P 

the  minimum  potential  energy  principle  requires  that 


3 L"P 

X; 


f tt  (x  , x_  . . . x ) ] = 0 i = 1,  2,  3,  . . . n 

p 1 2 n 


where  the  x are  the  unknown  discrete  parameters, 

i 
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From  the  expression  for  the  potential  energy  in  Appendix  C,  there  can 
be  formed  a set  of  matrices  by  rearranging  the  terms  as  follows  so  that  the 
resulting  matrices  are  symmetric  about  the  diagonal  of  each  matrix. 


where  each  vector  is  of  sice  5-2  “ and  each  matrix  is  of  total  size  5-2  x 5* 


which  is  divided  into  the  submatrices  as  indicated. 


The  first  matrix  (matrix  D)  contains  all  the  terms  associated  with  the 
potential  energy  of  the  sandwich  panel.  The  second  matrix  (E)  contains  all  of 
the  terms  associated  with  the  work  that  is  applied  to  the  sandwich  panel  from 
Appendix  C.  The  only  difference  is  that  the  applied  loads  are  now  written  as 
fo  llows : 


(N  _»N  »N  >N  „)=*•  (n  >n  R'n  *n  »n  r) 

xo  x-ts  xy  yo  y B xo  xB  xy  yo  yB 

where  now  the  submatrix  contains  the  relative  magnitudes  of  the  applied  loads, 

thus  being  able  to  factor  the  parameter  X out  of  the  matrix. 


To  apply  the  principle  of  minimum  potential  energy,  it  is  required  now 
to  take  the  partial  derivative  of  the  potential  energy  against  the  vector  X which 
contains  the  unknown  discrete  parameters.  When  this  is  done  the  resulting 
expression  will  be  as  follows: 

a* 

= O = DX  + X EX  (11) 
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11 

12 
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O 

42  Zx422 

2 2 
42  x i 
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2 X 4 2“ 

t2xr 

Now  the  expression  of  the  potential  energy  is  in  the  required  form  to 
apply  the  principle  of  minimum  potential  energy.  It  is  now  just  a matter  of 
solving  the  system  for  the  unknown  parameters: 

A , B , C , H , K , and  \ . 
mn  mn  mn  mn  mn 

This  is  easily  accomplished  and  will  be  shown  in  the  next  section. 

3.  3 ASSEMBLY  OF  EQUATIONS  FOR  SOLUTION 

To  easily  solve  the  system  of  equations  in  expression  (11),  the 
expression  must  first  be  rewritten  as  follows: 

DX  = - X EX  (12) 


To  solve  this  system,  one  can  multiply  the  submatrices  in  block 
form  as  follows  to  result  in  the  expressions: 


above.  This  is  accomplished  in  the  following  manner.  Rearranging  the  first 
expression  , it  can  be  rewritten  as  follows: 


* 


tv  'V 


V 'V 


and  now  solving  for  to  obtain 


'V  = - (V1  'V  'V 


Now  Equation  (14)  can  be  rewritten  as  follows: 


(14) 


(Xjl  = - [p]  [x^ 


(15) 


where  the  matrix  P is  the  multiplication  of  fD  ] * 

2 2 11 
a size  of  4^  x 1 . Therefore,  it  follows  that: 


f°12] 


and  therefore  has 


[P] 


-1 


tV  tD! 


(16) 


for  which  an  equation  solver  can  be  used  to  determine  the  entries  in  the 
matrix  P for  which  the  solution  is  needed.  Now  the  solution  for  the  vector 
needed  in  Equation  (15)  can  be  substituted  into  the  second  expression  in 
Equation  (13)  to  obtain: 

tv  1 -fPl)  [V  + fV  'V  = -x  f°]  'V- 

The  left  side  of  this  expression  can  be  rewritten  as  follows: 


<-tDu]  (Pl 


[V  1 


'V' 


The  matrix  in  the  above  expression  is  the  final  stiffness  matrix  which  is 
as  follows: 
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Therefore,  the  second  expression  in  Equation  (13)  can  be  written  as 


follows: 


where  all  the  entries  in  the  final  stiffness  matrix  K and  the  final  geometric 
stiffness  matrix  G result  from  the  equation  for  the  potential  energy  in 

Appendix  C. 

Finally,  the  critical  loads  and  the  corresponding  mode  shapes  of  the 
panel  can  be  obtained  by  solving  the  eigenvalue  problem  represented  by 
Equation  (17). 


SECTION  4 


SUMMARY  OF  RESULTS 

The  procedure  which  was  outlined  in  Section  3 has  been  implemented 
in  a computer  program.  Results  have  been  obtained  for  several  examples 
of  individual  critical  loads  from  other  articles  in  order  to  verify  the 
program's  accuracy.  Also,  critical  loads  have  been  calculated  for  various 
cases  of  combined  loads  and  these  results  have  been  compared  with  the 
existing  interactive  formulas. 

4.  1 EXAMPLES 

The  results  of  typical  calculations  performed  by  the  computer  program 
(FSCPAN)  to  test  its  validity  are  presented  in  Table  1.  The  results  are 
varied  basically  due  to  the  type  of  analysis  for  clamped  sandwich  panels 
performed  by  the  references  stated.  Each  reference  takes  a slightly  varied 
approach  to  the  analysis. 

For  these  examples,  agreement  with  the  previously  published  results 
was  good  in  the  case  of  shear  loads,  but  was  much  lower  in  the  case  of 
compressive  loads.  The  critical  loads  calculated  for  compression  are  much 
lower  than  those  cited  because  the  modes  given  do  not  uncouple.  Therefore, 
it  is  not  possible  to  solve  the  equations  in  closed  form  as  in  the  simple 
support  case  in  Reference  8. 

Convergence  of  these  examples  is  limited  to  eight  terms  where  £ = 8 
in  Section  3 (Equation  6)  of  this  report.  Using  eight  terms,  the  amount  of 
computer  space  needed  is  approximately  one- half  the  available  computer 
capacity  on  the  CYBER  74  at  Wright- Patter  son  Air  Force  Base.  If  needed, 
more  terms  could  be  used,  but  the  amount  of  space  needed  increases 
drastically  as  more  terms  are  added  to  the  expression  (Equation  6). 

As  for  pure  edgewise  bending,  there  were  no  previous  results  found 
to  verify  the  computer  program.  The  only  results  available  are  those  of 
the  interaction  formulas  which  will  be  examined  more  closely  later. 
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VERIFICATION  RESULTS 


4.  2 LEVY  ANALYSIS 


Critical  loads  calculations  for  typical  isotropic  panels  subjected  to 
compressive  loads  are  presented  in  Figure  4.  Samuel  Levy  (Reference  1) 
has  determined  a relation  between  the  critical  stress  ratios  and  the  length 
to  width  ratio  for  isotropic  flat  panels  subjected  to  an  edgewise  compressive 
load.  His  analysis  starts  with  the  case  of  the  simple  support  panel  in  which 
the  displacements  of  the  panel  are  assumed  to  be  zero;  then  by  revising  this 
analysis  he  is  able  to  adjust  the  coefficients  so  that  the  slope  of  the  clamped 
edges  is  also  zero.  His  analysis  leads  to  an  infinite  series  solution  for 
which  convergence  is  extremely  rapid.  In  Figure  4 his  results  are  plotted 
as  indicated  for  the  critical  stress  ratio 

(c  ) b2h 

x cr 


where 


12(1  - ^ 2) 

To  verify  the  analysis  of  this  report,  three  panels  were  used  to  obtain 
results  similar  to  that  of  Levy.  These  panels  are  presented  in  Table  2. 

The  panels  presented  were  vised  in  this  analysis,  computer  program 
(FSCPAN)  and  MIL-HDBK-23A  analysis  for  clamped  panels.  As  can  be  seen 
in  Figure  4 the  results  of  FSCPAN  arc  higher  than  the  Levy  analysis  due  to 
the  fact  that  the  FSCPAN  results  are  not  an  infinite  series  solution,  but 
are  lower  than  that  of  the  Handbook  analysis  for  the  length  to  width  ratios 
of  the  three  panels. 

For  the  Handbook  the  results  are  exact  solutions  for  the  three  panels, 
but  the  FSCPAN  results  are  for  eight  terms  of  the  infinite  series  needed  to 
solve  the  system.  With  more  terms,  the  results  will  be  much  lower  than 
those  presented  and  will  eventually  become  much  closer  to  the  Levy  results. 
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HANDBOOK 


Figure  4.  Levy  Critical  Stress  Ratio 


TABLE  2.  LEVY  PLATES 


r 


4.3  INTERACTION  FORMULAS 

A number  of  interaction  formulas  have  been  proposed  in  order  to 
facilitate  the  analysis  of  the  stability  for  cases  of  combined  loadings  for 
sandwich  panels  as  given  in  Reference  6.  The  purpose  of  the  interaction 
formulas  is  to  account  for  the  combined  effects  on  the  panel  of  various 
applied  loads  by  the  comparison  of  the  intensity  of  each  load  to  the  critical 
value  of  the  same  load  acting  alone  on  the  panel.  Each  ratio  is  weighted 


by  means  of  an  exponent  associated  with  the  particular  load  type  being 
examined.  For  example 


e e 

R 1 + R 2 = const. 

1 2 

where 

N. 

R i = N"  ’ 

icr 

The  e are  constant  exoonents,  and  the  index  L is  associated  with  the  type  of 
i 

loading  being  examined. 

Interaction  formulas  are  attractive  in  that  any  number  of  load  cases 
can  be  analyzed  once  the  individual  critical  loads  for  the  panel  are  determined. 
The  problem  is  that  these  formulas  represent  approximations  to  an  exact 
analysis  of  the  stability  under  the  combined  loadings.  Interaction  formulas 
for  the  cases  of  edgewise  compression  with  edgewise  shear,  edgewise 
compression  with  edgewise  bending,  and  edgewise  shear  with  edgewise 
bending  have  been  evaluated  for  the  sandwich  panels  in  Table  3 during  this 
analysis. 

4.  3.  1 Edgewise  Shear  and  Compression 

The  interaction  formula  for  axial  compression  and  edgewise 
shear  is  of  the  following  form: 
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TABLE  3.  PHYSICAL  DATA  FOR  COMBINED-LOAD  EXAMPLES 
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N 

x 

N 

xcr 


N 

/ N \ 

— L + 

r1  = >• 

N 

\ N / 

ycr 

xycr  / 

For  simplicity  of  this  analysis  and  for  the  purpose  of 
graphing,  the  above  formula  has  been  modified  to  the  following: 


2 


xcr  xycr 


The  left-hand  side  of  Equation  (18)  has  been  evaluated  for 
several  load  conditions.  The  resulting  values  of  the  equation  have  been 
graphed  and  are  presented  in  Figure  5 for  the  two  panels.  The  formula  is 
reasonably  accurate  for  most  of  the  cases  considered.  For  Panel  1 the 
estimate  (Equation  18  ) is  a nonconservative  estimate  in  that  the  true  values 
for  the  panel  are  less  than  the  estimate  below. 

4.3.2  Edgewise  Bending  and  Compression 

The  interaction  formula  for  combined  edgewise  bending  and 
axial  compression  loads  is  given  as  follows: 


N 

x_ 

N 

xcr 


N 


3/2 


xB 


1. 


N 


xBcr 


(19) 


Here  is  a pure-bending  load  as  shown  in  Figure  1.  The  left-hand  side 

of  the  formula  has  been  evaluated  for  several  load  conditions.  The  resulting 
values  of  the  equation  have  been  graphed,  and  are  presented  in  Figure  6 for 
the  two  panels.  The  formula  is  reasonably  accurate  for  most  of  the  cases, 
but  has  the  greater  discrepancy  for  Panel  1.  For  both  panels,  the  estimate 
(Equation  19)  is  a conservative  estimate  for  the  true  values  Qf  the  edgewise 
compression  and  edgewise  bending  for  the  panel. 


2f. 


CLAMPED  EDGES 


sion  - Shear  Interaction 


Compression 


4.  3.  3 Edgewise  Bending  and  Shear 

The  interaction  formula  for  the  evaluation  of  critical  com- 
bined edgewise  bending  and  shear  loads  is  as  follows: 
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xB 
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xBcr 


xy 


xyc  r 


= 1, 


(20) 


Values  of  the  left-hand  side  ol  the  formula  have  been  evaluated  for  the  two 
panels.  These  values  have  been  graphed  and  are  shown  in  Figure  7.  The 
expression  is  relatively  accurate  for  tin  two  panels.  For  both  panels  the 
estimate  (Equation  20)  is  a nonconscr\ ative  estimate  of  the  true  values 
for  the  edgewise  bending  and  shear  for  the  panels.  For  Panel  1 the  estimate 
is  almost  the  true  values  for  the  panel  in  that  the  graph  for  Panel  1 almost 
lies  on  top  of  the  formula.  For  Panel  2 the  estimate  is  much  more  non- 
conservative, in  that  the  true  values  for  the  panel  lie  much  farther  below 
the  estimate  for  the  values. 

4.4  OTHER  EXAMPLES 

Out  of  curiosity,  the  analysis  led  to  the  three  panels  listed  in  Table  4. 
For  each  of  the  three  panels  several  values  of  the  length  to  width  ratio  were 
considered.  Results  for  edgewise  compression  only  were  computed  for 
the  three  panels  with  the  computer  program  (FSCPAN)  and  the  Handbook. 

As  can  be  seen  in  Figures  8,  9,  and  10,  the  critical  edgewise  compression 
load  obtained  for  each  of  the  three  panels  is  lower  by  the  analysis  presented 
in  this  report  (FSCPAN)  than  by  the  analysis  presented  in  the  Handbook 
(Reference  6). 
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TABLE  4.  SANDWICH  PANELS 
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Figure  9.  Handbook  - FSCPAN  Comparison, 


75001- 


Handbook  - FSCPAN  Comparison. 


4.5  DISCUSSION 


In  this  section,  various  calculations  concerning  the  stability  of  a 
clamped  sandwich  panel  under  combined  modes  of  loading  and  single  loading 
have  been  presented  along  with  conclusions  concerning  the  effectiveness  of 
typical  interaction  formulas.  These  results  indicate  the  versatility  of  the 
analysis  procedure  described  herein,  and  illustrate  the  character  of  some 
of  the  effects  due  to  combined  edge  loadings  to  clamped  panels. 

The  results  given  in  Paragraph  4,  Z indicate  that  the  infinite  series 
solution  to  the  analysis  agrees  with  that  of  Levy  for  the  type  of  analysis  he 
has  performed.  Also,  for  the  eight-term  solution  used  for  this  analysis,  the 
results  for  an  edgewise  compressive  load  applied  to  an  isotropic  panel  are 
much  lower  than  those  results  as  presented  in  the  analysis  in  Reference  6. 
The  results  given  in  Paragraph  4.4  also  indicate  that  the  analysis  presented 
herein  for  an  edgewise  compressive  load  yields  lower  values  than  those 
values  as  presented  in  Reference  6. 

Results  given  for  the  various  interaction  formulas  indicate  that  the 
formula  proposed  for  edgewise  compression  and  edgewise  bending 
(liquation  19  ) may  be  of  some  use  for  design  purposes,  since  this  formula 
yields  rather  conservative  results  for  all  cases  considered  herein.  For  the 
case  of  combined  edgewise  compression  and  shear  loads  (Equation  18)  the 
formula  is  conservative  for  a rectangular  panel,  but  is  nonconservative 
for  a square  panel.  As  for  predictability,  the  formula  will  probably  depend 
on  the  type  of  panel  being  considered  and  the  length-width  ratio  being 
considered.  As  for  the  edgewise  bending  and  shear  interaction  formula 
(Equation  ZO  ),  the  results  obtained  will  be  nonconservative  for  the  design 
of  a sandwich  panel.  The  best  results  will  be  obtained  for  a square  sandwich 
panel  as  shown  in  Figure  7.  As  for  a rectangular  panel  the  formula  is  very 
nonconservativc  and  may  possibly  lead  to  a rather  large  error  in  the  true 
values  when  predicting  the  combined  load  values. 
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SECTION  5 


SUMMARY  AND  CONCLUSIONS 

An  analysis  has  been  presented  for  the  prediction  of  the  instability 
of  a flat,  clamped,  rectangular  sandwich  panel  loaded  by  combinations  of 
biaxial  edgewise  compression,  biaxial  edgewise  bending,  and  edgewise 
shear  loads.  The  model  has  been  formulated  in  terms  of  the  total  potential 
energy  of  the  panel,  which  is  composed  of  internal  elastic  strain  energy  and 
the  potentials  of  the  combined  applied  loads.  The  Ritz  method  is  used  to 
transform  the  energy  function  in  terms  of  continuous  field  variables  into  a 
quadratic  function  of  discrete  parameters.  Transforming  this  function  into 
a matrix  form  and  applying  the  principle  of  minimum  potential  energy 
results  in  a generalized  discrete  eigenvalue  problem.  Solution  to  this 
eigenvalue  problem  is  done  by  standard  methods.  The  finite  number  of 
eigenvalues  obtained  in  the  solution  approximate  the  lowest  value  of  the 
infinite  number  of  critical  loads  to  find  the  critical  buckling  load,  while 
the  eigenvectors  approximate  the  buckled  mode  shapes  that  the  panel  will 
assume  corresponding  to  the  given  eigenvalue. 

The  accuracy  of  this  analysis  has  been  demonstrated  by  the  comparison 
with  results  available  in  the  literature  for  a rather  limited  number  of  example 
Also,  a limited  study  has  been  made  to  look  at  the  applicability  of  several 
typical  interactive  formulas  which  are  commonly  used  in  the  design  of 
sandwich  panels.  From  the  present  analysis,  comparisons  of  the  interactive 
formulas  indicates  that  some  of  these  formulas  yield  nonconservative  results 
for  certain  combinations  of  loadings  and  geometries,  and  should  definitely 
be  used  with  extreme  caution. 

A computer  program  has  been  developed  to  implement  the  combined- 
loads  instability  analysis  presented  herein.  The  program  is  operational 
only  by  batch  mode  processing.  For  the  program  to  operate  interactively, 
the  number  of  terms  in  the  expression  (Appendix  A)  would  have  to  be  so 


r 


drastically  reduced  that  the  results  would  not  be  very  useful  to  obtain  the 

true  value  as  can  be  approximated  with  the  batch- mode  with  eight  terms. 

A detailed  description  of  the  computer  program  and  its  usage  are  contained 

. 9 

in  a companion  report. 

In  conclusion,  the  sandwich  panel  instability  analysis  described 
herein  represents  a useful  and  somewhat  accurate  tool  for  the  design  of 
light-weight,  high-performance  structural  components  having  clamped 
edges  which  is  the  case  in  many  physical  structures.  In  addition,  the 
method  of  analysis  provides  a suitable  starting  point  from  which  many  other 
general  types  of  sandwich  construction  can  be  considered. 
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This  Appendix  presents  the  expansion  of  the  Equations  (5  and  6) 
and  the  symmetry  needed  for  the  expansion  of  Equation  (2). 
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APPENDIX  C 

POTENTIAL  ENERGY  IN  QUADRATIC  FORM 

This  Appendix  presents  the  final  expression  for  the  potential 
energy  of  the  sandwich  panel.  This  includes  Equation  (8)  into  the  Equation 
(2)  and  Appendix  B where 
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